The r-mode instability: Analytical solution with gravitational radiation reaction 



Oscar J. C. Dias* 
Department of Physics, University of Waterloo, 
Waterloo, Ontario N2L 3G1, Canada 
and 

Perimeter Institute for Theoretical Physics, 31 Caroline St. N., 
Waterloo, Ontario N2L 2Y5, Canada 



in 
o 
o 



(N 

m 
> 
m 

o 

o 
in 
o 

■ 

Oh; 

6 ■ 



Paulo M. Sat 

Departamento de Fisica and Centro Multidisciplinar de Astrofisica - CENTRA, 
F.C.T., Universidade do Algarve, Campus de Gambelas, 8005-139 Faro, Portugal 

(Dated: February 2, 2008) 

Analytical r-mode solutions are investigated within the linearized theory in the case of a slowly 
rotating, Newtonian, barotropic, non-magnetized, perfect-fluid star in which the gravitational radi- 
ation (GR) reaction force is present. For the GR reaction term we use the 3.5 post-Newtonian order 
expansion of the GR force, in order to include the contribution of the current quadrupole moment. 
We find the explicit expression for the r-mode velocity perturbations and we conclude that they are 
sinusoidal with the same frequency as the well-known GR force- free linear r-mode solution, and that 
the GR force drives the r-modes unstable with a growth timescale that agrees with the expression 
first found by Lindblom, Owen and Morsink. We also show that the amplitude of these velocity 
perturbations is corrected, relatively to the GR force-free case, by a term of order fi®, where Q, is 
the angular velocity of the star. 

PACS numbers: 04.40.Dg, 95.30.Lz, 97.10.Sj, 97.10.Kc 



I. INTRODUCTION 

Stars are not rigid systems; they naturally oscillate 
and the non-radial pulsations generate gravitational ra- 
diation (GR) that removes energy and angular momen- 
tum from the star. In non-rotating stars this process is 
always dissipative and, inevitably, the oscillation of the 
star dies off. However, if the star rotates this is not nec- 
essarily the case: under certain conditions the amplitude 
of the pulsation mode grows and a GR instability sets 
in. In this paper we are interested on a special class of 
GR-driven instabilities called Chandrasekhar-Friedman- 
Schutz (CFS) instabilities, after the pioneering works of 
Refs. [1, 2]. The generating mechanism of these CFS in- 
stabilities is well understood. Indeed, consider a mode 
that counter-rotates, i.e., whose propagation direction is 
opposite to the star's rotation in the co-rotating frame of 
the star. This mode has then negative angular momen- 
tum. When analyzing the properties of the mode in the 
inertial frame, one sees that the mode is dragged forward 
by the stellar rotation and can become prograde, i.e., it 
can rotate in the same sense of the star. Then, the gravi- 
tational waves emitted due to the mode oscillation carry 
away positive angular momentum from the star. Since 
positive angular momentum is extracted from a mode 
with negative angular momentum, the amplitude of the 
mode increases and an instability appears. So, although 
the GR emission lowers the inertial frame energy, under 
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the above conditions it induces an increase in the mode 
energy as measured in the co-rotating frame. These two 
events can occur simultaneously because the energy in 
the co-rotating frame E and the energy in the inertial 
frame E\a are related hy E = E^^ — VlJ , where Q, is the 
angular velocity of the star and J is the angular momen- 
tum of the mode. Therefore, E can grow if both E-^ and 
fiJ decrease [2]. The condition for the appearance of a 
CFS instability can be stated as follows. Given a mode 
with an angular pattern speed a in the inertial frame that 
propagates in a star that rotates with angular velocity fi, 
a CFS instability sets in if < fi < il. In general, this 
means that for a given mode with fixed a the instability 
will appear only if the star is rapidly rotating with an 
angular velocity greater than the critical value f2crit — c. 
This requirement of a large stellar rotation, together with 
the unavoidable viscosity contributions that damp the in- 
stability, has diminished the astrophysical interest in the 
CFS instabilities [3]. 

This pessimistic scenario changed when the properties 
of a special class of stellar pulsation modes, called r- 
modes, were investigated in detail. The r-modes, which 
were first studied [4-7] more than twenty years ago, are 
pulsation modes of rotating stars that have the Coriolis 
force as their restoring force. These modes induce per- 
turbations mainly in the fluid's velocity and cause very 
small disturbances in the star's density. More precisely, 
the velocity perturbations are of order fi, while the den- 
sity perturbations are of order fi^, where VI is the star's 
angular velocity. The unique feature that clearly distin- 
guishes the r-modes from other stellar modes and that 
turn them potentially interesting for astrophysics is the 
fact that they are driven unstable by GR reaction in all 
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perfect-fluid rotating stars, no matter how slowly they 
rotate. This property was discovered by Andersson [8], 
and afterward confirmed more generally by Friedman and 
Morsink [9]. So, unlike other modes, r-modes are always 
retrograde in the star's co-rotating frame and prograde 
in the inertial frame, i.e., the sign of the r-mode fre- 
quencies is always opposite in the two frames. In other 
words, the CFS instability always occurs for the r-modes 
because the angular pattern speed of a r-modc with an- 
gular quantum number I in an inertial frame is given by 
cr = ^'7(V+i|^^ ^ [4] ' ^^'^ ^° satisfies the CFS condition, 
< a < f2, for any value of the stellar rotation and for 
/ > 2. 

Two questions are then naturally raised, namely, what 
are the timescales associated with the growth of the r- 
mode CFS instability, and are the several viscosity pro- 
cesses and damping mechanisms present in a star suf- 
ficient to suppress the GR instability? These questions 
were first addressed in Ref. [10], where, for the evaluation 
of the growth timescale t of the instability, the small den- 
sity perturbations induced by the r-mode were neglected 
and, in the framework of Newtonian hydrodynamics, the 
formalism developed in Ref. [11] was used to compute the 
GR reaction that induces the GR instability. GR cou- 
ples to r-modes primarily through the current multipoles 
rather than the usual mass multipoles (this result follows 
straightforwardly from an order count in the multipoles: 
the perturbations in the mass multipoles are of order fP, 
while the perturbations in the current multipoles are of 
order Q). In Ref. [10] it was found that the energy E of 
the r-mode perturbation in the co-rotating frame varies, 
to lowest order in O, according to i;^ = £^06"^*^^, where 
Eq is the initial energy. Moreover, for / = 2, the timescale 
tgr of the GR-driven instability is given by [10] 

where G is the Newton's constant, c is the velocity of 
light, and J = drpr^ with p and R being the density 
and the surface's radius of the unperturbed star, respec- 
tively. In parallel, a similar programme as the one of 
Ref. [10] was carried in Ref. [12], but with the impor- 
tant difference that, in order to include the small den- 
sity perturbations, the calculation was performed up to 
second order in the stellar rotation (fi^). The density 
perturbation contribution is important to correctly eval- 
uate the bulk viscosity effects, but has no influence on 
the GR timescale tgr [12]. The calculations performed 
in Refs. [10, 12] indicate that, for a wide range of rel- 
evant temperatures and angular velocities of the star, 
TGR ^ Ty, with the gravitational timescale tor be- 
ing several orders of magnitude smaller than the viscous 
timescale ry. Thus, for this range of temperatures and 
angular velocities of the star, viscosity cannot suppress 
the growth of the GR-driven instability of r-modes. 

The astrophysical implications of the GR-driven insta- 
bility in the r-modes are then very promising. Indeed, 



a long-standing unsolved problem of pulsar astrophysics 
is associated with the observed slow rotation rates of 
neutron stars. Neutron stars are believed to be formed 
through stellar gravitational collapse. The angular mo- 
mentum of the system is conserved during this process 
and so the natural expectation is that the final neutron 
star should have a rotation rate close to the Keplerian 
velocity, i.e., the maximum velocity above which matter 
starts escaping through the equatorial plane (see never- 
theless Ref. [13] for a different view). However, all the 
observed neutron stars present a rotation rate well below 
this maximum value. The r-modcs with their peculiar 
properties provide a possible explanation for these small 
rotation rates of young pulsars in supernova remnants. 
Indeed, as we just saw, in a newly born, hot, rapidly- 
rotating neutron star the GR reaction force dominates 
bulk and shear viscosities for enough time to allow most 
of the star's angular momentum to be radiated away as 
gravitational waves [10, 12]. As a result, the neutron 
star might spin down to just a small fraction of its initial 
angular velocity. The r-modes may also play an impor- 
tant astrophysical role in low- mass X-ray binaries. In- 
deed, GR emission due to the r-mode instability could 
balance the spin-up torque due to accretion of neutron 
stars in low-mass X-ray binaries, limiting in this way the 
maximum angular velocity of these pulsars to values con- 
sistent with observations [14-16]. GR released during 
the r-mode evolution seems also to be a very promis- 
ing source for the detection of gravitational waves, even 
though the perturbations in the density induced by r- 
modes are small and thus, in a first intuitive analysis, 
they should not be associated with strong GR emission. 
A first attempt to develop a time-dependent model of the 
evolution of the r-mode instability, including an approx- 
imate discussion on the nonlinear phase where most of 
the angular momentum of the star is radiated away, has 
shown that the GR emitted can probably be detected in 
enhanced versions of laser interferometer detectors [17]. 
Although recent results on the nonlinear saturation of 
the r-mode energy are not so optimistic, they still point 
to the possible detection of gravitational waves if the un- 
stable neutron stars are close enough [18]. 

A quite interesting feature that has emerged from re- 
cent investigations on r-modes is the presence of differen- 
tial rotation induced by the r-mode oscillation in a back- 
ground star that is initially uniformly rotating. That 
differential rotation drifts of kinematical nature could be 
induced by r-mode oscillations of the stellar fluid was 
first suggested in Refs. [19, 20], where an approximate 
analytical evaluation of nonlinear effects was performed. 
Differential rotation drifts induced by r-modes were also 
found in numerical simulations of nonlinear r-modes car- 
ried out both in general relativistic hydrodynamics [21], 
and in Newtonian hydrodynamics [22-24]. Differential 
rotation was also reported in a model of a thin spherical 
shell of a rotating incompressible fluid [25] . Recently, an 
analytical solution, representing differential rotation of 
r-modes that produce large scale drifts of fluid elements 
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along stellar latitudes, was found within the nonlinear 
Newtonian theory up to second order in the mode am- 
plitude and in the absence of GR reaction [26]. This 
work is exact to second order, improving the approxi- 
mate analytical computation performed in Refs. [19, 20]. 
As shown in Ref. [27], this differential rotation plays a 
relevant role in the nonlinear evolution of the r-mode in- 
stability. Indeed, in a newly born, hot, non-magnetized, 
rapidly rotating neutron star, r-modes saturates a few 
hundred seconds after the mode instability sets in; the 
saturation amplitude depends on the amount of differen- 
tial rotation at the time the instability becomes active 
and can take values much smaller than unity [27]. Now, 
the question that remains unsolved is whether GR reac- 
tion also induces differential rotation, making an extra 
contribution. One of the aims of the present paper is to 
initiate a programme that hopefully will allow to answer 
this question. 

In this paper, we will perform an analytical study 
of r-mode solutions within the linearized theory in the 
case of a slowly rotating, Newtonian, barotropic, non- 
magnetized, perfect-fluid star in which a gravitational ra- 
diation reaction force is present. The GR reaction contri- 
bution is assumed to be given by the 3.5 post-Newtonian 
order expansion that includes the current quadrupole mo- 
ment, which is the main responsible for the GR instabil- 
ity that sets in. The detailed formalism of the 3.5 post- 
Newtonian order expansion is given in Refs. [28, 29]. It is 
known that the most unstable mode is the I = m = 2 r- 
mode, where / and m are the angular momentum and the 
azimuthal numbers, respectively. Hence, we will focus 
our attention on the GR-driven instability of this mode. 
We will find an analytical expression for the r-mode ve- 
locity perturbations, which is physically consistent with 
the energy evolution studied in Refs. [10, 12]. 

The plan of the paper is the following. In Section II we 
write the full Newtonian hydrodynamic equations with 
the GR reaction force. We review the main properties 
of the GR reaction force and we give the general for- 
mulas needed to compute this reaction force. The ex- 
plicit computation of the GR reaction force in the 3.5 
post-Newtonian order expansion is then performed in Ap- 
pendix A. Section III is devoted to the linearized theory. 
The perturbed Newtonian hydrodynamic equations with 
the GR reaction force are written in Subsection III A. 
In Subsection III B we briefly review the well-known GR 
force-free linear r-mode solution that is used in Subsec- 
tion III C to compute the first-order Eulerian change in 
the GR force. In Subsection HID we finally solve the 
Newtonian hydrodynamic equations with the GR force. 
In Section IV our results are discussed and future direc- 
tions are pointed. In Appendix B, we identify explicitly 
the connection between some of our intermediary for- 
mulas for the GR force and those of previous numerical 
works on the subject. 



II. NEWTONIAN HYDRODYNAMIC 
EQUATIONS WITH GR REACTION FORCE 

The Newtonian hydrodynamic equations for a uni- 
formly rotating, barotropic, non-magnetized, perfect- 
fluid star in the presence of the gravitational radiation 
(GR) reaction force are the Euler, continuity, and Pois- 
son equations given, respectively, by 

dtV + {v-V)v = -p-^VP~V<i> + F'^^, (2) 

dtp + V ■{pv) = 0, (3) 
= AnGp, (4) 

where p is the density of the fluid, v is its velocity, P is its 

pressure, $ is the Newtonian potential, and F'^^ is the 
GR reaction force per unit mass (or, simply, GR reaction 
force). The general expression for the GR reaction force 
i^*^^, that includes contributions from the time-varying 
mass multipole moments and from the time-varying cur- 
rent multipole moments, has been found and progres- 
sively discussed in detail up to the 3.5 post-Newtonian 
order. A review of this path and the 3.5 post-Newtonian 
expansion of the GR force can be found in Refs. [28, 29]. 
In this paper we are interested on the contribution of 
the GR reaction force to the r-mode oscillations. Con- 
trary to what happens in other mode oscillations, the 
r-mode instability is predominantly excited by the time- 
varying current multipole moments rather than by the 
usual time- varying mass multipole moments [10] (see also 
a detailed discussion in Ref. [29]). Now, in order to 
include the contribution of the current multipole mo- 
ments, the post-Newtonian expansion of the GR reac- 
tion force must be done at least up to the 3.5 order 
[28, 29]. This post-Newtonian order includes the con- 
tribution of the current quadrupole moment, which is 
the dominant current multipole moment. Since the main 
contribution for the r-mode instability comes from the 
current quadrupole moment, we will henceforth neglect 
all the terms coming from the time-varying mass multi- 
pole moments, i.e., we will not consider terms that appear 
in post-Newtonian orders below 3.5 neither the terms of 
the 3.5 post-Newtonian order that have their origin in 
the time- varying mass multipole moments. These terms 
are clearly identifled in Refs. [28, 29]. So, in the so-called 
Blanchct's gauge, the contribution of the current multi- 
pole moments to the 3.5 post-Newtonian order GR force, 
henceforth labelled as for simplicity, is ^ 

F^^ = -dJ+vx{V xP), (5) 

where /3 is a vector whose components are given by 



^ See equations (12) and (17) of Ref. [29]. Our notation is slightly 
different from the one of Refs. [28, 29]. The correspondence is 
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with I = 1, 2, 3 and similarly for the other Latin indices. 
In the previous expressions, Sij{t) is the time- varying 
current quadrupole tensor, 

Stj{t) = J (fixekq(iXj)XkpVg, (7) 

Cijk is the Levi-Civita tensor, is the Cartesian coordi- 
nate of the point at which the tensor is evaluated, and 

sl?it)^^s,M (8) 

The explicit computation of the GR reaction force is per- 
formed in Appendix A. 

It is interesting to note that /? can be appropriately 

identified as a gravitational vector potential [25]. This 
designation comes from the clear equivalence between 
the GR force per unit mass exerted on a moving fluid 
(5) and the Lorentz force per unit charge exerted by an 
electromagnetic field on a moving charged fluid, F^°'^ = 
—dtA + v X (V X A), where A is the electromagnetic vec- 
tor potential. In this analogy, the designation of gravita- 
tional vector potential for /? follows straightforwardly if 
one identifies (3 = A. 

III. i?-MODES IN THE LINEARIZED THEORY 
WITH THE GR REACTION CONTRIBUTION 

In previous investigations [10, 12], the hydrodynamics 
equations with the GR force were used to obtain an ex- 
pression for the time evolution of the physical energy of 
the r-modc perturbation, dE/dt, from which the gravita- 
tional radiation and viscous timescales were determined. 
In this section, we explicitly solve the hydrodynamics 
equations with the GR force, obtaining an expression for 
the first-order Eulerian change in the velocity of the r- 
modes S^^^v. It is shown that these velocity perturbations 
are sinusoidal with the same frequency as the well-known 
GR force-free linear r-mode solution; that the GR force 
drives the r-modes unstable with a growth timescale that 
agrees with the expression found in Refs. [10, 12]; and 
that the amplitude of these velocity perturbations is cor- 
rected, relatively to the GR force-free case, by a term of 
order fl^. 



where we have defined 

^(i)C7=^+5(i)$. (12) 
P 

In the above equations, 

V = QrsmOe^ (13) 

is the velocity of the unperturbed star and p its mass den- 
sity, while S^^^Q denotes the first-order Eulerian change 
in a quantity Q. 

B. Linear r-mode solution without the GR 
reaction contribution 

In our approach, we use the linear r-mode solution 

without the GR contribution to generate the GR reac- 
tion force. As already mentioned, the GR force drives 
r-modes unstable and the most unstable mode is the 
I = m = 2 r-mode. Wc will therefore concentrate our 
attention in this mode and in this subsection we briefly 
present this well-known solution. 

In the absence of the GR reaction force, S^^^'F'^^- — 0, 
and to lowest order in a and in O, where a is the dimen- 
sionloss r-modc amplitude, the perturbed equations (9)— 
(11) allow r-mode solutions with velocity perturbations 
given by 5^^^v — aflR{r / RYYif e^"^ , where R is radius of 
the unperturbed star. In spherical coordinates and for 
1 = 2, the velocity perturbations are explicitly given by 

^(^V = 0, (14a) 

5i%, = -L.[I^r^sinee'^^^+^'\ (14b) 

4 \ TT R 

<5(% = l./I^r^ sine cos ee^(^*+'^*\ (14c) 
and from Eq. (9) one obtains 

^(i)f/= i,/^^r3sin2ecosee*(2^+-*), (15) 

6 V TT it 

where oj = —40/3 is the I = 2 r-mode frequency in the 
inertial frame [4]. 



A. Linearized hydrodynamic equations 

The hydrodynamic equations (2) (4) for a uni- 
formly rotating, Newtonian, barotropic, non-magnetized, 
perfect-fluid star with the GR reaction force can be lin- 
earized yielding, in the inertial frame, 

= -V,5(i)[/ + 5Wi?™, (9) 
dtS^^'ip + 1! ■ V(5(iV + V • (p(5(iV) = 0, (10) 
y^W^ = 47rG5(iV, (11) 



C. Explicit computation of the first-order Eulerian 
change in the GR force 

To find the linear r-mode solution of the perturbed 
hydrodynamic equations (9)-(ll) with the GR reaction 
force wc have to evaluate the first-order Eulerian change 
in the GR force, i.e., we have to compute (5(^'F^^. To 
do so we expand both the velocity vector, v = v + 
and the current multipole tensors, Sij = Sij + S^^^Sij, 
where Sij describes the multipole tensor of the unper- 
turbed star and (J^^'S'jj is the first-order Eulerian change 
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in the multipolc tensor. The first-order Eulerian ehange 
in the GR force S^^^F'^^ contains then terms of the type 

The multipolc tensors of the unperturbed star are eas- 
ily computed from Eqs. (A3), with the use of Eq. (13), 
yielding 



0: 



t,j = x,y,z. 



(16) 



i.e., they all vanish. In practice, this result allows us 
to rewrite the first-order Eulerian change in the GR 
force simply as a sum of terms of the type Vk5^^'>sf^ and 

Now, to compute the first-order Eulerian change in the 
multipolc tensors S^^'^Sij we use an approach in which we 
assume that the linear r-mode solution with = 

acts as a source for the current multipole tensor. In other 
words, and taking S^^^Sxx as an example, we assume that 



/ 



dV sin 6 cos ( 



S'^^^piyo sin (f) + v^ cos 6 cos > 



+ p[S^^^V0 sin (j) + S^^^v^ cos 6 cos ( 



.(17) 



where v is given byEq. (13) and (5( ^">v is given by Eqs. (14) 
with uj being now an arbitrary parameter to be deter- 
mined. Since S^^^p is of order afl^, S^^^v is of order aft, 
and V is of order 17, in Eq. (17) there are terms of order 
aft, namely, pS^^'^Vi, and terms of order a^l^, namely, 
ViS^^^p. Since our interest is to find the r-mode solution 
induced by the GR force to lowest order in a and in 
CI, we will neglect in Eq. (17) the contribution coming 
from the terms ViS^^^p and we will work only with the 
dominant terms p5^^\. This discussion applies similarly 
to the other current multipole tensors, whose first-order 
Eulerian change can be straightforwardly written from 
Eqs. (A3). 

Under the above conditions, Eqs. (14) are inserted into 
Eq. (17) yielding 



(18) 



where uj is an arbitrary parameter, 0{aCl'^) denotes terms 
of order afl^ or more, and J is defined as 



J= drpr^. 



(19) 



Similarly, it is straightforward to find that the first-order 
Eulerian change in the other multipole tensors satisfy the 
relations 

yy — "-^xx"} 
SWSx, = = d^'K. = 0. (20) 



From Eqs. (18) and (20) it is then clear that the n^^ time 
derivative of the multipole tensor perturbation is given 

by 



(21) 



In spite of the large number of terms that appear in 
the definition of the GR force [see Eqs. (A4) in the Ap- 
pendix A], Eqs. (18), (20), and (21) allow a considerable 
simplification of the first-order Eulerian change in the 
GR force. Indeed, to lowest order in O, the first-order 
Eulerian change in the GR force is given simply by 

j(i)i^GR ^ -3inian^iv''r^smHcos9e'^^^+^'\ (22a) 



~ iK^aOo;^ (uj + 30 sin^ 6) sin 06'^^'''+"'^ , 
R 



(22b) 



SWpGR ^ -K^aQu'^r^ sine cosee'^^'^+'^'\ (22c) 

where the constant k sets the strength of the GR reaction 
force and is defined as 



16 /ttG 



K = 



45 V 5 c7 



(23) 



D. 



Linear r-mode solution with the GR reaction 
contribution 



We have now all that is needed to find the linear r- 
mode solution that solves, to lowest order in CI, the lin- 
earized fluid equations (9)-(ll) with the GR force. 

It is well established that the GR drives the r-modes 
unstable. This means that the GR reaction induces a 
small complex imaginary part in the allowed r-mode fre- 
quency, i.e., one has 



(24) 



where ujq = Re[aj] is the frequency of the r-mode and 
w = Im[ix;] < is the small imaginary part that is related 
to the growth timescale of the instability of the mode. 

Assuming that w/ujo ^ 1, an assumption that will 
be checked at the end, the perturbed GR force (22) can 
finally be written as 



^WpGK ^ aClojl (wo + 30sin2 9) sin 6* 



GR 



J 



(25a) 



(25b) 



-K-aClcolr^ smecosee-'^*e'^^'^+'^°*l (25c) 
R 
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The perturbed Euler equation (9), in the inertial 

frame, is given by 

(26a) 



GR 

I 



(26b) 

+ ^a^5(^)c/ = (5(^)i^«^, (26c) 
rsin^ ^ 

with defined by Eqs. (25). In order to solve 

Eqs. (26), we try a solution that has the same structure 
as the GR force-free r-mode solution given by Eqs. (14) 
and (15). More specifically, we try an Ansatz of the type^ 



0, 



(27a) 



5 aQ. 



^sin6'e-'^*e*(2'^+'^ot) 



W0 



4 V TT i? 

+ 5(%(i,r,e,.^), (27b) 
1 . /l!!^^2 sin 6 cos 6 e— *e*(2^+-ot) 

4 V TT i? 

+ (5(%(t,r,0,(/.), (27c) 



and 



1 /5c^ 
6V7r R 



sin2 6'cos6'e-'^*e'(2'^+'"°*^ 



(28) 



where J^^V and J'-^'C/ are functions of i, r, 9 and </> to be 
determined, and w and wq are unknown constants to be 
fixed. 

As already mentioned above [see Eq. (1)], it is known 
from previous investigations [10, 12] on the time evolu- 
tion of the physical energy of the r-mode that w = 1/tgr 
is proportional to ffi. This result follows straightfor- 
wardly from an inspection of Eqs. (26) with S^^^F^^ and 
S^^^v defined by Eqs. (25) and (27), respectively. In- 
deed, the time derivative of S'^^^v on the left-hand side 
of Eqs. (26) yields terms proportional to Qzu, while the 
right-hand side of these equations contains, to leading or- 
der in O, terms proportional to O^. Similarly, the r-mode 



frequency uq could also have a correction (relatively to 
the GR force- free value coq = — 4f2/3) proportional to Q^. 
Therefore, we write w and coo in the form 



and 



Wo 



(29) 



(30) 



where a and b are (by now arbitrary) real dimensionless 
constants and we have defined 



,10 



J 



(31) 



Note that for a = 1 the growth timescale of the GR- 
induced r-mode instability, tgr = l/w, coincides with 
the one obtained in Rcfs. [10, 12]. For 6 = 0, the fre- 
quency coo coincides with that of the GR force-free linear 
r-mode solution. 

For such an Ansatz the perturbed Euler equation re- 
duces to the following system of differential equations for 
<5(i){;^ and S'^^'^lf, 



—20 sin 

= ijan'^r^ sin^ 6 cos ee-^*e'(2^+'^°*) , 



(32a) 



= ij 



-{a + bi + 2) -sin^i 



ailJr'^ sm0 



X g-I37tgi(2(^-|-Wot) 



(32b) 



{dt + nd^)d'^'^\ + 2ncoseS^^^ve 
27 



r smt 



= --^{a + bi + 2) aQ'^r'^ sin cos 6le-'^*e'(2'^+'^''*' . 

(32c) 

The right-hand side of the above system induces a so- 
lution of the form 

= i7aO*^rV(^) e-'^*e'(2'^+'"°*\ (33a) 
= -jan^r^g{e)e-^'e'^'^'^+^°'\ (33b) 



and 



<5(i)i7 = ljan''r^h{e)e-^*e'^^'^+^'>'\ 



(34) 



At a first glance we might expect that the GR r-mode solu- 
tion should have a non- vanishing radial velocity component since 
the first-order Eulerian change in the GR force (25) has a non- 
vanishing radial component. However, this is not necessarily the 
result that we might expect, as the experience acquired with the 
GR force-free r-mode indicates. Indeed, in this case there is no 
GR force, but the Coriolis force, which is the restoring force, 
is present. Now, although the Coriolis force has a radial com- 
ponent, the fact is that, to leading order, the force-free r-mode 
velocity perturbations (14) do not have a radial component. 



where f{9), g{9) and h{9) are functions of 9 determined 
by the system of equations 



2 sin Og{0) + ih{6) = i sin^ 9 cos 9, 
2f{9) - 6cos9g{9) - ideh{9) 



-3i sin 9 



- {a + bi + 2) - sin^ I 



3ism9cos9f{9) -ism9g{9) - h{9) 
= -\{a + bi + 2) s].T?9cos9. 



(35a) 
(35b) 

(35c) 
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Prom the system of equations (35) it is straightforward to 
obtain a hnear first-order differential equation for h{9), 
namely, 

dh(e) l-3cos2 6», 
dO smycosy 

where the source term q{6) is given by 



q[9) = -{a + bi- l)sin( 



(37) 



This equation has the solution 
h{e) = ^sin^6icos6i 

+ -(a + 6i- l)sin^6icos6'ln(tan6'), (38) 

where A is a (complex) constant determined by initial 
data. In Eq. (38), the first term of the right-hand side 

corresponds to the general solution of the homogeneous 
equation [Eq. (36) with q{0) = 0], and the second term on 
the right-hand side corresponds to the particular solution 
of Eq. (36) with the source term q{0). From Eqs. (35) 
and (38) is then straightforward to obtain: 

f{e) = -^(A- l)sin6i-|- ^(a + 6i- l)sin6» 

- — (a-|-6i- I)sin6iln(tan6i), (39) 
g{e) = -^(A- 1) sin 6* cos 6* 

- — {a + bi- 1) sine cose \n{ta.n 6). (40) 

The above solution to the perturbed Euler equation 

must also satisfy the continuity equation (10), which at 
order afl^ is simply given by 



V- (5(1)^=0. 



(41) 



Inserting S'^^^v from Eq. (33), with f[9) and g{d) given 
by Eqs. (39) and (40), into the continuity equation (41) 
we obtain: 



[smef{0)]-2g{e)=O 



(42) 



or 



4i sill 2z 
- — (a + bi - 1) + —{a + bi - l)sin6icos6i = 0,(43) 

9 cosO 9 

implying that a = 1 and 6 = 0. 

Therefore, we arrive at the conclusion that the r-mode 
velocity perturbations are sinusoidal with the same fre- 
quency as the well-known GR force-free linear r-mode so- 
lution, and that the GR force drives the r-modes unstable 
with a growth timescale that agrees with the expression 
found in Refs. [10, 12]. 

Inserting solutions (33) and (34) [with f{0), g{0) and 
h{9) given by Eqs. (38)-(40), a = 1 and 6 = 0] into 



Eqs. (27) and (28) we finally get that the r-mode solution 
to the linearized fluid equations with the GR force is 
given by 



(44a) 



tan 



r sine 



X g-^i'tg'(2<^-|-wot) 



(44b) 



V4, 



-an 



5 1 



2 V 7ri? 



ij{A - i)n^ 



r smycos( 



(44c) 



and 



-an^ 



1 /5 1 



2 V TTi? 



+ ijAn^ 



sin^ cos ( 



X g-'3'*g»(2<^-|-wot) 



with zu and wq given by 



and 



LOO 



8 hr 



jRn'^ 



4Q 



(45) 



(46) 



(47) 



The velocity perturbations given by Eqs. (44) have a 
piece similar to the GR force-free solution, the differ- 
ence being the factor e~^* responsible for the exponen- 
tial growth of the r-modc amplitude due to the presence 
of a GR reaction force, and another piece proportional to 
a'y{A — l)n^, where A is a constant fixed by the choice 
of initial data. If this constant is chosen to be of order 
unity, A ~ 0{1), then the second term in the square 
brackets of Eqs. (44) and (45) is much smaller than the 
first term. Indeed, let us assume that the mass density 
p and the pressure P of the perfect-fluid star are related 
by the usual polytropic equation of state P = kp^, with 
k such that M = IAMq and R = 12.53 km. For such a 

choice of the equation of state, J = 1.12 x 10^^ kg m^ and 

the Kcplcrian angular velocity at which the star starts 
shedding mass through the equator is Ox ~ •^y/nUp = 
5612 sec^^, where p = 3.38 x 10^^ kg/m"^ is the aver- 
age mass density of the star and G is the gravitational 
constant. Then, in Eqs. (44) [or similarly in Eq. (45)] 
the e xpression inside the square brackets can be written 
as v/5/(47ri?2)[l + l.Qg x 10-^i{A - l)(0/0;f)5]. For 
A ~ 0(1) and n ^ n^, the second term in the square 
bracket can be neglected. 

Note that in the GR force- free limit, the constant k 
that sets the strengh of the GR reaction force goes to 
zero. Therefore, in this limit 7 defined in Eq. (31) goes 
to zero and the S^^^'v perturbations vanish as they should, 
since they are proportional to 7. 
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IV. DISCUSSION OF THE RESULTS 

Analytical r-mode solutions were investigated, within 
the linearized theory, in the case of a slowly rotat- 
ing, Newtonian, barotropic, non-magnetized, perfect- 
fluid star in which a gravitational radiation (GR) reac- 
tion force is present. For the GR reaction term wc have 
used the 3.5 post-Newtonian order expansion of the GR 
force, in order to include the contribution of the cur- 
rent quadrupole moment. The evaluation of the GR 
force perturbation requires the knowledge of the multi- 
pole moment perturbation and its time derivatives. To 
compute this Eulerian change in the multipole moment 
we have used an approach in which we assumed that the 
GR force-free linear r-mode solution acts as a source for 
the current multipole moment. We have carried our anal- 
ysis of the perturbed fluid equations only to lowest order 
in fi. This order is enough to evaluate the GR reaction 
contribution. A higher order analysis would include per- 
turbations in the density which would be important if we 
were interested in the accurate evaluation of the viscosity 
effects. 

We have found the analytical expression of the r-mode 
velocity perturbations that solves the Newtonian hydro- 
dynamic equations with the GR reaction force. Prom this 
solution, given by Eqs. (44) and (45), four important fea- 
tures of the r-mode instability driven by the GR reaction 
can be read: 

(i) The velocity perturbations S^^'^v are proportional to 
^ii24>+ujot) ^ with ivo = -in/3. Thus, they have the same 
sinusoidal behavior and the same frequency wq as the 
solution of the GR force- free linear velocity perturbations 
given by Eqs. (14). 

(ii) The amplitude of the velocity perturbations is pro- 
portional to exp{— tut}. Since < 0, the GR force in- 
duces then an exponential growth in the r-mode ampli- 
tude. The e-folding growth timescale tgr = l/w agrees 
with the GR timescale (1) first found in Ref. [10]. Thus, 
these velocity perturbations are consistent with the ex- 
pression found in Ref. [10] for the time evolution of the 
energy of the r-mode perturbation, dE/dt. 

(iii) The velocity perturbations 5^^\] contain a piece pro- 
portional to 7(^4 — l)a.n^ , where A is an arbitrary con- 
stant fixed by the choice of initial data. If we choose this 
constant A to be of order unity, then this part of the 
solution could be neglected. 

(iv) The parameter k, defined in Eq. (23), sets the 
strength of the GR reaction force. Thus, the GR force- 
free limit is obtained when we set k equal to zero. In 
this limit, w and 7 defined, respectively, by Eqs. (29) 
and (31) go to zero. Then, from Eqs. (44) and (45) we 
recover the GR force- free linear r-mode solution given by 
Eqs. (14) and (15). 

Our results arc strictly valid only for the slow-rotation 
regime since our solution holds only to leading order. 
However, there are good indications that the evolution 
of r-modes in rapidly-spinning stars does not differ sig- 
nificantly from the evolution in slowly-rotating stars [30] . 



Hence, it is reasonable to expect that at least the quali- 
tative features that we found are also valid in the rapid- 
rotation regime. 

A natural extension of the present work is to try to find 
an analytical r-mode solution of the nonlinear hydrody- 
namic equations with the GR reaction force. The nonlin- 
ear GR force- free case was analyzed in Ref. [26] , where an 
analytical r-mode solution in the nonlinear theory up to 
second order in the mode amplitude was found. This so- 
lution represents differential rotation that produces large 
scale drifts of fluid elements along stellar latitudes. It 
contains two separate pieces, oiic^ induced by first-order 
quantities and another determined by the choice of ini- 
tial data. Since these two pieces cannot cancel each 
other, differential rotation is an unavoidable nonlinear 
kinematic feature of r-modes. The analysis of Ref. [26] 
can be extended to the case in which the GR reaction 
force is present. The main purpose of this extension is to 
find if the GR reaction force provides an extra source of 
difii'erential rotation. This work is in progress [31]. 
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APPENDIX A: EXPLICIT COMPUTATION OF 
THE GR REACTION FORCE 

In the cartesian basis the explicit components of the 
gravitational vector potential (6) are 




The use of the usual relations between cartesian (a;, y, z) 
and spherical (r, 9, (f)) coordinates yields the components 
of the gravitational vector potential (3 in spherical coor- 



dinates, 

Pr = 0, 

(3e = -^r^ cotiO (cos 045] 



45c 



sm 



+ sin( 



Syy = J dVpr^ sin 9 sin (p {v$ cos (j) — v^ cos 6 sin cf)) , 
(A2a) (A3b) 
Szz = J dVpr^v^ sin 6 COS 6, (A3c) 

cos(2<^)5lfl + ^ sin(2(/.) _ 45] j I ^ (A2b) S,y = ^ J dVpr^ [ve sin ^ cos(2<^) - v^, sm{20) sin(2<^)] , 

(A3d) 

Sxz ~ ~2 [ ^^P'''^ 6sm(f> + cos(2^) cos (f)] , 

ism{2e) [Spj + Sg - 2SS + cc.s(2^,) (sfj - Sg)] (^3^, 

+ 2cos(2.)3in,49+si„(2.)3M2,)sg}, (A2c) 5„ = 1 / .V^..^". - . c„. , - c„(2.) .„ * (A3t, 

where the current multipole tensor in spherical coordi- 
nates is given explicitly by 



8G 



45c 



+ 



Sxx = — J dVpr^ sin 9 cos (j) {vq sin (fy + v^ cos 6 cos </>) 



Inserting Eqs. (A2) into Eq. (5) yields the spherical 
(A3a) components of the GR reaction force, 



7.GR 



_16G 
~15^' 

8G 
+ 15^™^^ 



rve |cos0 (cos05| - sin</.5S) + sin0 ^ sin(2<^) (-5^ + + cos(20)5i:; 



I sin(2^) sin(2<^)5[fl + ^ sin(2^) [s (4^1 + 45j) + cos(2</.) (4^1 - ^^j) 
+2cos(26i) (cos(/)4^J +sinfl;>4l) 



(A4a) 



16G 
45^' 



AG 



r < cos f 



(cos ct^Sf, - sin05[fj) + sin0 cos{2ct>)S^^l + i sin(20) (4^^ _ 



+ 



45c7 



jrvA'i (4'i + 4?) (1 + 3cos(20)) - 6sin2 6 [cos(2<^) - ^Ifj) + 2sin(2<^)5^1 



16G 
+ ^ cos 



-12sin(26') (cos(/)41 + sin^4^) 

(cos<^45 - sin ,^4^1) + sine \\ sin(2.^) (-4^1 + 4^j) + cos(2</,)4^ 



(A4b) 



8G 



8G 



{29) cos<pSfl + i sin(2e) [s (46I + 46J) + cos(20) [sfl - 4^^ 
+ 2 cos(26l) sin (^5|fJ + sin(26l) sm{2(l))Sfl | 
^g^,, sin(20) sin(20)45] + \ sm{29) [s (4^1 + 45J) + cos(20) (451 - Sf^ 



AG 
"45c7' 



+2cos(2e) (cos.^5^J + sin 68]'. 
rve^^ (^i + 41) (1 + 3 cos(20)) - 6 sin^ 9 [cos(2</.) (4^1 - S^^l) + 2 sin(2</.)45l 
-12sin(2e) (cos</)5^J +sin</,4^l) |, 



(A4c) 
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where we have made use of the relation Szz = ^{Sxx + 
Syy) [see Eqs. (A3)] to remove Szz from the final expres- 
sion. 

As a check to our expressions for the GR reaction 
terms, we apply them to the spherical shell system dis- 
cussed in detail in Ref. [25] . The authors have computed 
the instability timescale driven by the GR reaction and 
concluded that it is given by Eq. (1) with J = pR^, and 
p and R being the surface density and the radius of the 
spherical shell. Now, if we introduce the r-mode velocity 
given by Eqs. (1) and (7) of Ref. [25] into our Eqs. (A2) 
and (A3), we obtain Eqs. (B14) and (55) of Ref. [25] that 
define the timescale of the instability on the shell. 



APPENDIX B: CONNECTION TO PREVIOUS 
WORKS 



In previous works, in which the effects of the GR force 

on the r-mode cvohition were studied numerically, the fi- 
nal form of the GR force was always written as a function 
of the multipole moments Jim, instead of the quadrupole 
tensors Sij [22-24]. To compare our expressions with 
those of Refs. [22-24] we establish the connection be- 
tween the two approaches in this Appendix. 

The current multipole moments are defined as [11] 



Jim= / dVpr'v-Yi 



(Bl) 



where I and m are the angular momentum and the az- 

imuthal numbers, respectively, Yi^* is a spherical har- 
monic vector of the magnetic type defined by [11] 

YiZ' = (-1)™^-™ (B2) 

and 

Y,^ = [l{l + l)]-'/'rx{VYim). (B3) 
The usual spherical harmonic functions Yim are given by 



Yim = i-iy 



2l + l{l-m)l 



-,1/2 



47r {l + m)\ 



e'""l'Pr {cose), (B4) 



where 6 and (j) are the angular spherical coordinates and 
P™ (cos 9) are the associated Legendre functions defined 
by 

1 fjl+m 

Pr(cose) = -7- sin"* 9- ^7^(cos2 9 - iV. (B5) 

Like in our case, the authors of Refs. [22-24] are inter- 
ested only in the most unstable mode, the / = 2 r-mode. 
This mode excites the J2m current multipole moments 



which in a spherical coordinate system are written as 



J22 
J21 



-^"''r'^ sin 9 {v^ cos 9 + ive) , (B6a) 
\\ - f dVpe-''^r^ [v^ cos(26') + ive cos 9] , 

4 V TT J 



J: 



20 



/ '^^/'^'^^sin(20). 



(B6b) 
(B6c) 



From Eqs. (A3) and (B6) it is straightforward to ver- 
ify the following correspondence between the current 
quadrupole tensor Sij and the current multipole mo- 
ments J2m [23]: 



Syy Sx 



2iSxy = 4 



J2 



(B7a) 
(B7b) 
(B7c) 

The most unstable r-mode is in fact the I = m = 2 

mode which excites J22 but not J21 a-nd J2o. Hence, in 
Refs. [22-24], the GR reaction force is computed assum- 
ing that, in Eq. (B7), J21 — and J20 — 0, i.e., 



SxZ — iSyz — 2\l —J2 



^xx ~t~ ^vv — ^zz — ^1 



-'VV 



Sxx "I" 2iS'_ 



xy 



|J22, 



Sxz ~ iSyz — 0, 

^xx ~l~ ^yy ~ Szz — 0* 

It is also assumed there that 

4^1 ~ (iw)"J22, 



(B8a) 

(B8b) 
(B8c) 



(B9) 



and it is found that the GR force in the inertial frame is 
given by 



^GR 



3K{ive — cos 9)r sin 9e^'^'^ J22 , 
-mr^ sin 06^*^4? 



(BlOa) 



22 



+ 3K{-iVr + V4, sm9)r sin6ie^''^4f , (BlOb) 
FpR ~ Kr^ sin 9 cos 9e^"^J^^ 



+ 3K{vr COS 9-ve sin 9)r sin 9e^"^J^^ . (BlOc) 

Now, let us apply the perturbative approach used in 
our paper to check if our results match with Eq. (BIO). If 
we linearize the multipole moments, Jim = Jim Jim, 
where Jim describes the multipole moment of the unper- 
turbed star and 6^^^Jim is the first-order Eulerian change 
in the multipole moment, it is straightforward to show, 
using Eqs. (B6) and (13), that 



Jim=0; Z = 2;m = 0,l,2, 



(Bll) 
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and the use of Eqs. (B6) and (14) yields 



aO, — e" 
R 



together with 



(B12a) 

(B12b) 
(B12c) 



^(1)4^1 = («C^)",5(1)J22. (B13) 
In fact, the value of k in Refs. [22-24] is two times ours (compare 
our Eq. (23) with Eqs. (A.8) and (A. 9) of Ref. [23]). Their extra 
factor of 2 appears to be a typographical mistake. Indeed, with 
our value of K we find an exact agreement with the instability 
timescale of the r-mode in a spherical shell [see discussion after 



So, the right-hand side of the first-order Eulcrian change 
of Eqs. (B7b) and (B7c) is zero in agreement with tlic 
assumptions (B8b) and (B8c) taken in Refs. [22-24], and 
the derivative of 5^^^J22 satisfies (B13) in agreement with 
the approximation (B9) assumed in Refs. [22-24]. More- 
over, if we insert the value of S^^^J22 given in Eq. (B12a) 
into the first-order Eulerian change in the GR force we 
obtain an expression that agrees with our result (22) '. 



Eq. (A4)] and we also get for the Newtonian star a GR growth 
timescale that agrees with the one found in Refs. [10, 12]. If we 
use the value of k defined in Refs. [22-24] we obtain a timescale 
that is twice the correct one in both geometries. 
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